fl and f2 + production in central Pb+Pb collisions at 40 and 158A GeV 
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Results are presented on Q production in central Pb+Pb collisions at 40 and 158A GeV beam 
energy. For the first time in heavy ion reactions, rapidity distributions and total yields were measured 
for the sum fl~+Cl + at 40 A GeV and for f2~ and f2 + separately at 158A GeV. The yields are 
strongly underpredicted by the string-hadronic UrQMD model but agree better with predictions 
from hadron gas models. 
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The measurement of multi-strange particles is of par- 
ticular interest in heavy ion collisions at ultra-rclativistic 
energies. One important aspect is the observation that 
the inverse slope parameter T of the f2 mt spectrum is 
significantly smaller than expected from the linear mass 
dependence of T naively implied by the presence of ra- 
dial flow. This led to the hypothesis that multi-strange 
hyperons are not affected by the pressure generated by 
the hadronic matter in later stages of the reaction [2|. 
Originally, the increase of the production of multi-strange 
particles as compared to elementary hadron-hadron colli- 
sions was suggested as a signature of quark-gluon plasma 



formation [3j. However, existing experimental data on 
EE and A production at lower beam energies 0, ex- 
hibit a much stronger enhancement than observed at top 
SPS energies. Generally, it is found that the abundances 
of strange particles are close to those calculated in sta- 
tistical models assuming the creation of an equilibrated 
hadron gas In a hadronic environment, as expected 
at lower beam energies, this equilibration is generally 
difficult to achieve. At larger energy densities, when 
the hadronic system might be close to the QGP phase 
boundary, multi-particle fusion processes could lead to 
fast equilibration 0. However, there exists no dynamic 
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explanation in a hadronic scenario at lower energy densi- 
ties. The present measurement of f2 at 40^4 GeV provides 
an important test for these models. Recent results on the 
energy dependence of the ratio (K + )/(ir) 0, || indicate 
a sharp maximum of relative strangeness production at 
a beam energy of 30A GeV. This observation can be in- 
terpreted as a signal for the onset of deconfinement 01 
and might be reflected in the energy dependence of multi- 
strange particle production. 

The data were taken with the NA49 large acceptance 
hadron spectrometer at the CERN SPS. With this detec- 
tor, tracking is performed by four large-volume TPCs. 
A measurement of the specific energy loss dE/dx pro- 
vides particle identification at forward rapidities. Time- 
of-flight detectors improve the particle identification at 
mid-rapidity. Centrality selection is based on a measure- 
ment of the energy deposited in a forward calorimeter by 
the projectile spectators. A detailed description of the 
apparatus can be found in [llj. 

We present in this paper an analysis of two samples 
of central Pb+Pb events taken at beam energies of 40 
and 158A GeV in the years 1999 and 2000, respectively. 
About 5.8 x 10 5 events were recorded at 40^4 GeV with a 
centrality selection of 7.2% of the total inelastic cross sec- 
tion corresponding, on average, to (N w ) — 349 wounded 
nucleons At 158 A GeV, 2.8 x 10 6 events were taken 
at 23.5% centrality corresponding to (N w ) — 262. 

The f2 were identified in the decay channel — > AK, 
A — » p7r (68% branching fraction). To reconstruct the fV 
(f2 + ), the A (A) candidates were selected in an invariant- 
mass window of 1.101-1.131 GeV/c 2 and combined with 
all negatively (positively) charged tracks in the event. 
The same procedure as in the S analysis of |13| was used 
to identify the secondary vertex of the ft decay. 

To reduce the combinatorial background several cuts 
were applied. Identification of the (anti-)protons by 
dE/dx in the TPCs reduced the contribution from fake 
A (A). The measured dE/dx was required to be within 
3.5 standard deviations from the predicted Bethe-Bloch 
value. Likewise an enriched kaon sample was extracted 
from the charged tracks. A further background reduction 
was achieved by requiring a minimal distance of 25 cm 
in the beam direction between the target and the d de- 
cay vertex position. The candidates were extrapolated 
back to the target plane to obtain the transverse coordi- 
nates b x (magnetic bending plane) and b y of the impact 
point with respect to the primary interaction vertex. To 
reject non- vertex candidates, cuts of \b x \ < 0.5 cm and 
\b y \ < 0.25 cm were applied. Kaons from the primary 
vertex were excluded by imposing a cut of \b y \ > 1.0 cm 
on the kaon tracks. In addition, \b y \ > 0.4 cm was re- 
quired for the A candidates at AO A GeV. With these cuts 
an acceptable separation of signal and background was 
achieved for fl transverse momenta above 0.9 GeV/c. 

In Fig. the invariant-mass distributions of the fl~ 
and candidates are shown for p t > 0.9 GeV/c and 
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FIG. 1: The invariant-mass distributions of and Q, + can- 
didates. Left: summed distribution of AK - and AK + pairs 
at 40A GeV. Middle: AK" pairs at 158A GeV. Right: AK+ 
pairs at 158A GeV. The full curves represent a fit to signal 
and background described in the text. The dashed curves 
show the background contribution. 



—0.5 < y < 0.5. Note that the available statistics at 
40 A GeV is not sufficient to separately analyze the 
and + [l4|. Clear signals are observed at the fl mass 
of mo = 1672.5 MeV/c 2 0] with a resolution of 5 and 
4 MeV/c 2 at 40 and 158A GeV, respectively. 

The spectra were fitted to the sum of a polynomial 
background and a signal distribution, determined from 
the simulation described below. The raw Q yield is ob- 
tained by subtracting the fitted background in a mass 
window of ±7 MeV/c 2 around the nominal mass. 

Detailed simulations were made to correct the yields 
for geometrical acceptance and losses in the reconstruc- 
tion. For this purpose, a sample of Q was generated in 
the full phase space accessible to the experiment. The 
Geant 3.21 package 0] was used to track the generated 
f2 and their decay particles through a detailed description 
of the NA49 detector geometry. NA49 specific software 
was used to simulate the TPC response taking into ac- 
count all known detector effects. The simulated signals 
were added to those of real events and subjected to the 
same reconstruction procedure as the experimental data. 
The acceptance and efficiency were calculated in bins of 
Pt and y as the fraction of the generated f2 which tra- 
verse the detector, survive the reconstruction and pass 
the analysis cuts. This fraction amounts in total to 0.4- 
0.5% (0.2-0.5%) for 40 (158)^ GeV data, depending on 
p t and y. The geometric acceptance is of the order of 
20%, which in turn is reduced to 0.6-1.2% by the cuts 
that suppress the combinatorial background. A further 
reduction of the reconstruction efficiency at 158^4 GeV 
by 30-60% is due to the high track density. 

The statistical error is given by the quadratic sum of 
three contributions: the signal, the background and the 
efficiency correction, the latter being smaller at 40 than 
at 158A GeV. The systematic uncertainties are domi- 
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nated by the background subtraction method and by im- 
perfections in the simulation. By varying the analysis 
strategy and the cuts applied, a systematic error of 10% 
is estimated in the transverse- mass region (ra t — mo) > 
0.3 GeV. At lower m t this error is about 25% 

In Fig. |21 the transverse-mass spectra of the S7 are 
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FIG. 2: The transverse-mass spectra of S7 and Sl + at mid- 
rapidity. Left: SI - + Cl + at 40A GeV (triangles) and fi~ at 
158A GeV (circles). Right: H+ at 158A GeV. The errors 
shown are statistical only. The dashed curve shows the expo- 
nential fit described in the text. The full, dotted, and dash- 
dotted curves show a model including transverse expansion. 

shown integrated over the range ±0.5 (40^4 GeV) and 
±1 (158 A GeV) around mid-rapidity. The transverse- 
mass spectra were fitted to an exponential function 

dN / m t \ 

«exp -— (1) 



m t dm t dy 



in the range (mt — mo) > 0.2 GeV. The results are plotted 
in Fig. |21 (dashed curves) and the inverse slope param- 
eters T are listed in Table |U No significant difference 



TABLE I: The inverse slope parameter T (MeV), the width 
a of the rapidity distribution, the mid-rapidity yield dN/dy 
and the total yield (N) of SI production at 40 and 158A GeV. 
The first error is statistical and the second systematic. 



40A GeV 


158A GeV 


158A GeV 






si+ 


T 218±39±39 


267±26±10 


259±35±18 


a 0.6±0.1±0.1 


1.2±0.4±0.2 


1.0±0.4±0.2 


dN/dy 0.10±0.02±0.02 


0.14±0.03±0.01 


0.07±0.02±0.01 


(N) 0.14±0.03±0.04 


0.43±0.09±0.03 


0.19±0.04±0.02 



can be observed in the shape of the SI - and Sl + spectra 
at 158A GeV. The inverse slope parameters are close to 
the values obtained by the WA97 and NA57 collabora- 
tions 0, ^3 ■ The inverse slope parameter at 40 A GeV is 
somewhat lower but, within errors, compatible with the 
158 A GeV result. 



To investigate whether the SI decouples earlier from the 
fireball than lighter hadrons, we use a hydrodynamical 
model which assumes a transversely expanding emission 
source 18]. The parameters of this model are the freeze- 
out temperature Tf and the transverse flow velocity /3 S 
at the surface. Assuming a linear radial velocity profile 
fit(r) = f3 s r/R, which is motivated by hydrodynamical 
calculations, the mtspectrum can be computed from 



dN 



m t dm t dy 



, 'pt sinh p\ ( mt cosh p\ 
rdrm t I [ — — \KA — I, (2) 



where R is the radius of the source and p = tanh _1 /3 t is 
the boost angle. The full curve (A) in Fig. [21 shows the 
result of a calculation with Tf = 90 MeV and an average 
flow velocity (/?t) = 0.5. These parameters were obtained 
from a simultaneous fit of the model to the mtspectra of 
K + , K~, p. p, ^ A, and A, all measured by NA49 at 
158A GeV [IlEllllIii- The dotted curve (B) is calcu- 
lated with T f = 170 MeV and (A) =0.2, obtained from 
a fit to J/ip and tp' spectra j21|. The disagreement of 
curve (A) and the agreement of curve (B) with the data 
suggest that in this version of the model, the freeze-out 
conditions of the $2 are similar to those of the J/-0 or tp' 
but are different from those of the lighter hadrons. The 
use of a constant expansion velocity, on the other hand, 
results in a fair agreement of the measured hadron spec- 
tra, including the S and the SI, with Tf = 127 MeV and 
(f3 t ) = 0.5 .22], as shown by curve (C) for the Omega. 
It predicts, however, a significant decrease of the spec- 
trum for heavy hadrons for m t — mo — » 0. This dip is 
already quite pronounced for the Omega but not sug- 
gested by the measurements [Isj. Thus, a radius inde- 
pendent transverse expansion velocity may be too crude 
an approximation of the velocity profile. 

The parameterizations of Eqs. and J2J were used to 
extrapolate the SI yields into the unmeasured regions of 
m t . Assuming that the shape of the mt distribution does 
not depend on rapidity, extrapolation factors of 2.3 (2.2) 
at 40 (158)A GeV were obtained from fits to the summed 
£!~ and Sl + data. A systematic uncertainty of 6% is due 
to the choice of parameterization. 

The extrapolated SI yields at 40 and 158^4 GeV are 
shown in Fig. |3| as a function of rapidity. All spectra can 
be described by a Gaussian with zero mean and a width 
a obtained from a fit to the data, see TableHl The widths 
of the O - and Sl + spectra at 158A GeV are compatible 
but are both significantly larger than the width measured 
at 40^4 GeV. Also given in Table |U are the mid-rapidity 
yields dN/dy (-0.5 < y < 0.5) and the total yields (N) 
obtained by extrapolating the rapidity spectra into the 
unmeasured region using the Gaussian fits. The mid- 
rapidity yields are slightly below the values given by the 
NA57 collaboration [2j| , however, they agree within sta- 
tistical errors, if the difference in the centrality selection 
at 158^4 GeV is taken into account. 
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FIG. 3: The rapidity dependence of f2 production in central 
Pb+Pb collisions. Left: Q~ at 40 A GeV (triangles) and 

fi _ at 158A GeV (circles). Right: Q + at 15&4 GeV. The 
errors shown are statistical only. The open symbols show the 
measured points (full symbols) reflected around mid-rapidity. 
The curves correspond to Gaussian fits to the data. 



In the following we denote by (O) the sum of the to- 
tal and f2 + yields and by (it) the total charged pion 
yields from [j| , multiplied by a factor 1.5. The pion yields 
at 158A GeV were scaled by the ratio of the numbers of 
wounded nucleons to account for the difference in the 
centrality selection of the pion and the f2 measurement 
(note the centrality selection at AOA GeV of 7% and at 
158A GeV of 23.5% most central events.). Figure^Jshows 
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FIG. 4: The ratio {Q)/(n) (see text) versus the center-of-mass 
energy. Statistical errors are shown as lines, while the brackets 
denote the systematic errors. The dashed curve shows the 
prediction from the hadronic string model UrQMD |24jl and 
the gray box that of RQMD |2fij . A hadron gas model without 
strangeness suppression |2^ | is shown by the full curve. The 
open squares represent the fits from | 2q | including strangeness 
under-saturation. 



2.3 including the color rope mechanism [25j. 

On the other hand, the data are close to the predic- 
tions of statistical hadron gas models which use a grand 
canonical ensemble. In these models, the chemical freeze- 
out temperature and the baryonic chemical potential are 
fitted to the yields of other measured hadrons. The 
hadron gas model of [2^| (labeled B in Fig. 0J) introduces 
in addition a strangeness undersaturation parameter 7 S 
in the fits, which have been performed at each energy 
separately. The fit results at 30^4 GeV reflect the sharp 
maximum of the K + /7r + -ratio observed around this en- 
ergy. The present measurement at 40^4 GeV seems to 
favor this model, compared to that of j2]j (labeled A 
in Fig. which does not allow for strangeness under- 
saturation (7 S = 1 for all energies). The data point at 
158A GeV, however, does not discriminate between the 
two models. Nevertheless, the observation that f2 pro- 
duction is compatible with phase-space undersaturation 
at higher SPS energies (> 30^4 GeV), would be in line 
with a similar behavior of the kaon excitation function in 
the same energy regime 0. 

In summary, NA49 has performed a measurement of O 
production in central Pb+Pb reactions over a wide region 
of phase space. At a beam energy of 158A GeV the avail- 
able statistics allowed to separately analyze fl~ and f2 + . 
The shapes of the transverse-mass spectra at this energy 
reveal no difference between and £! + and are in agree- 
ment with previous results by WA97 and NA57. In a hy- 
drodynamically inspired model with radially increasing 
velocity profile, the data favor a low transverse expansion 
velocity and high freeze-out temperature. The rapidity 
spectra of the f2, which have not been measured before 
in heavy ion reactions, are compatible with a Gaussian 
shape. The widths for f2~ and f2 + appear to be similar. 
The yields are strongly under-predicted by the string- 
hadronic UrQMD model. The data agree better with 
predictions from hadron gas models. 
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the ratio (£l)/{ir) as function of the center-of-mass energy. 
The ratio tends to increase with energy. It is clearly un- 
derpredicted by the UrQMD string-hadronic model [24| 
as shown by the dashed curve in Fig. A better descrip- 
tion of the 158 A GeV data is provided by RQMD version 
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